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The detection of gravitational waves from compact binaries relies on a computationally bur- 
densome processing of gravitational-wave detector data. The parameter space of compact-binary- 
coalescence gravitational waves is large and optimal detection strategies often require nearly re- 
dundant calculations. Previously, it has been shown that singular value decomposition of search 
filters removes redundancy. Here we will demonstrate the use of singular value decomposition for a 
composite detection statistic. This can greatly improve the prospects for a computationally feasible 
rapid detection scheme across a large compact binary parameter space. 
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I. INTRODUCTION 

Ground-based laser-interferometric gravitational-wave 
detectors have demonstrated sensitivity to gravitational- 
wave strain at the level of 10 -22 or better over a band 
from 80 — 1000 Hz [1—3]. Later this decade advanced 
detectors will surpass the present sensitivity by a fac- 
tor of ~ 10 [4, 5]. One of the most promising sources 
of gravitational waves is the merger of two compact ob- 
jects [ ]. Advanced detectors are expected to detect on 
the order of 40 neutron star - neutron star merger events 
per year [7, 8] and a similar number of mergers involving 
black holes. The exact parameters of the signals, such 
as mass and spin of the component objects, will not be 
known ahead of time. Therefore, the optimal detection 
strategy must include the possibility of detecting signals 
with unknown parameters. Matched filtering has been 
employed to search this parameter space [9, 10]. The pa- 
rameter space is explored by choosing a discrete set of fil- 
ters that guarantees that all signals within the parameter 
space are found with a signal-to-noise ratio (SNR) greater 
than ~ 97% [11] of the maximum possible. Recently it 
has been shown that the use of singular value decomposi- 
tion can reduce the number of filters necessary to search 
the parameter space [12]. This paper extends that work 
to explore use of the singular value decomposition (SVD) 
filter outputs for detection without reconstructing the 
physical template waveforms. 

If one is interested only in knowing whether any of 
the signals are present, and not which one is present, 
then the problem of detection decouples from that of pa- 
rameter estimation. Wainstein and Zubakov (1962) [13] 
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describe the problem of detecting any of several signals 
without concern for parameter estimation as composite 
detection. This paper will explore the composite de- 
tection of compact binary signals using the techniques 
proposed in [12]. We find that the composite detection 
statistics explored produce a lower detection efficiency 
at a fixed false alarm rate than traditional approaches. 
However when combined hierarchically with traditional 
approaches, the combined approach can perform as well 
at low false alarm rate, but with reduced computational 
cost. 



II. COMPOSITE DETECTION OF COMPACT 
BINARY SIGNALS 

We will consider an unknown gravitational wave signal 
arising from a compact binary coalescence in the digitized 
gravitational-wave detector output as a vector of data 
points parameterized by an unknown amplitude A and a 
collection of unknown parameters 9 (e.g., the masses of 
the component bodies) that determine the shape of the 
signal. We will denote this signal as As(9). Both the 
amplitude and shape parameters are not known a priori. 
The relative frequency of parameters in the population 
is described by the probability density function p(A, 6). 
We will assume the joint distribution p(A, 0) is separa- 
ble, that is p(A, 8) = p(A)p(9). This is not generally 
true globally across the parameter space, but should be 
roughly true locally. We will denote the vector of dis- 
cretely time-sampled strain data as h. In addition to 
containing normally distributed noise, n, h will possibly 
contain the gravitational-wave signal, 

h = n + As(d). (1) 

Assuming the noise has unit variance and s(6) is nor- 
malized such that the inner-product of it with itself is 
unity, A is the expected value for the SNR of a signal 
after optimally filtering. 
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The optimal detection statistic is the marginalized like- 
lihood 



Using the change of basis described in [12], where 



A 



P(h|s) 
p(h|0)' 



(2) 



where p(h|s) is the probability of obtaining h given the 
presence of any signal, with the signal parameters A, 9 
integrated out 

p r+oo 

p(h\s)= / P (h\A,9)p(9)p(A)dAd9, (3) 

and p(h|0) is the probability of obtaining h in the absence 
of any signal. The marginalized likelihood increases as 
the probability of the data containing a signal increases. 
For white noise, 



p(h|0) oc exp ( ~^h T h 



(4) 



This white noise form is also valid for colored noise cases 
if one applies the linear whitening transformation to the 
data h. In the frequency domain, this transformation 
divides the data by the amplitude spectral density of the 
noise. Our goal is to find a computationally inexpensive 
approximation to this optimal detection statistic using 
the SVD-reduced filter set described in [12]. 



A. Expansion of the marginalized likelihood 

In this section we will consider an expansion of the 
marginalized likelihood and show that it gives rise to a 
detection statistic that exploits the SVD filter basis. 

For simplicity we will assume that the signal parame- 
ters 9 take on discrete values [14] so that we may index 
different signals as Si. We will also assume that the am- 
plitude, A, does not depend on the signal index i. The 
detector output now has the form 



h = n + Asi. 



(5) 



We don't know a priori which signal Si is present in the 
data. The marginalized likelihood [9] is 



A cx expL4h • sj], 



(6) 



where the integral in (3) is replaced by a sum over all pos- 
sible templates Sj. Expanding the exponential to second 
order yields 



l + Ah-s j + ^ 2 (h-Sj) 



(7) 



The first order term is oscillatory and contributes little 
to the sum. We thus ignore it and define the approximate 
likelihood as A' 



A' :=5>-Sj 



(8) 



(9) 



the approximate likelihood expression (8) simplifies to 

A' = ]TKh. Uk ) 2 , (io) 

k 

due to the properties of the orthonormal matrix vjk. (10) 
has no dependence on the original parameter index j. As 
was shown in [ I 2] , fewer filters are required using the u k 
basis, therefore the sum over k has fewer terms than the 
original sum over j. Unfortunately higher order terms in 
the expansion of (7) do not simplify as well. It is possible 
to derive detection statistics that perform better than 
(10) while still exploiting the reduced filter set; one such 
example is shown in the next section. 



Assuming the signal probability distribution is 
a multi-variate normal distribution 



In this section we explore approximating the signal 
probability p(h|s) in (3) as a multivariate normal distri- 
bution. This will result in a different detection statistic 
that still exploits the SVD filter basis. 

We begin by assuming the signal probability (3) has 
the form 



p(h|s)cxexp ^-ih T C- x h^ , 



(11) 



whose covariance matrix C of second-order moments 
completely defines it. The covariance matrix is 

Cij = cov(hi,hj) 

= ((m + As0)+){n 3 + Asj(0))) 
= (n inj ) + (A 2 ) (8i(§)8j(§)) 

= S i:i + (A 2 ) f Si {9) Sj {9)p(9) d0, (12) 
Jv 

where (A 2 } — J^°° A 2 p(A) dA and we have assumed the 
independence of the variables A, 9 and n^. 

Using (4) and (11), we compute the logarithm of the 
marginalized likelihood, A. This is monotonic in A and 
therefore is sufficient to use as a ranking statistic 



In A oc T 



h T I 



(13) 



The performance of (13) depends on the waveform pop- 
ulation; we will assess it for simple but realistic cases in 
Sec. III. 

We compute C by drawing N samples of 9 distributed 
according to p(9). We assume a signal population that 
is distributed according to our ability to distinguish sig- 
nals, which is approximated by the filter banks described 
in [11]. This definition of the signal population is ap- 
proximately uniform in the 9 distribution for sufficiently 
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small regions of parameter space. The integral in (12) is 
then approximated by a summation over these samples, 



a.. — a 

~ l J 



(A 2 ) N 



N 

fc=l 

2\ N 



N ^ 



: (0fc)s(0fc), 



(14) 



fc=i 



where are discrete values of 9 and s(9k) is a row vector 
of the time samples Sj(#fe). By defining the matrix Sik = 
Si{0k) we can simplify the notation to 



I 



<^> T 



N 



S A S. 



(15) 



This is precisely the arrangement of the signal matrix 
proposed in [12]. As in (9), we use the SVD to decompose 
the signal matrix S (here written in matrix notation but 
equivalent to (9)), 



S = VSU J 



(16) 



where V and U are unitary matrices and S is a diag- 
onal matrix of the singular values of S, referred to by 
components in the Sec. II A. U is the matrix of the 
orthonormal basis vectors described in the previous 
section. Then we note that 



i-cr 1 = i- i + ^s T 



N 



N 



US 2 U T 



= U I — 



N 



U J 



(17) 



where we have used UU T = U T U = I. Let us define a 
diagonal matrix J 



Jfcfe := 1 - 1 



(A 2 ) 
~2\T 



k 



,2 , JV_ 
T k T {A 2 } ' 



(18) 



With this definition the approximate detection statistic 
(13) can be written as 



T = h T UJU T h, 

which, in the notation of the previous section, is 

r 2 



r = E ^ /(A2) (h-u k ) 



(19) 



(20) 



It is worth noting the limits of this expression for low 
and high amplitude signals (small and large values of A) , 



I™ r = E^( h - u k) 2 > 



lim T = V(h-u k ) 2 . 



(21a) 



(21b) 



The A — > limit is exactly the same result as the sec- 
ond order approximation derived in the previous section. 
In this limit the most important filters (i.e., those with 
larger singular values) contribute more to the composite 
detection statistic. The A — > oo limit in (21b) reduces 
the composite detection statistic to the sum of squares 
of the orthogonal filter outputs, (h • Uj). This detection 
statistic is equivalent to the excess power statistic ob- 
tained in Anderson et. al for the detection of waveforms 
of known bandwidth and duration [15, equation (2.10)]. 
Here, instead of projecting the data onto a basis that 
spans a time-frequency tile, we project the data onto a 
basis that span the space of CBC waveforms. The essen- 
tial difference between (20) and the excess power statistic 
of Anderson et. al is that (20) folds in knowledge of the 
relative probability that the waveforms we seek match 
any of the basis vectors individually, whereas the target 
waveforms in [15] are assumed to match the basis vectors 
of the time-frequency tile with equal probability. 

Fig. 1 presents an example of the components of (18) 
for the small and large amplitude limits, and for an ampli- 
tude of 20. These were produced using the signal matrix 
described in Sec. Ill A. 



III. OPERATING CHARACTERISTICS OF THE 
PROPOSED COMPOSITE DETECTION 
STATISTICS 

In this section we explore the performance of the pro- 
posed detection statistic (20). We begin by establish- 
ing the framework with which we conduct simulations to 
produce Receiver Operator Characteristic (ROC) curves. 
These indicate the probability of detection versus the 
probability of false alarm. We then present some practi- 
cal scenarios in which to understand these results. 

We begin by assuming the detector data h has the form 
of (5) modified to allow an ambiguous phase 



h = n + ^4s i(0 ) + -Bs; 



i(7T/2)- ( 22 ) 

Under ideal situations the signal S; could be recovered by 
a matched filter with an SNR of \J A 2 + B 2 . According 
to [12] the signal Si can be decomposed into orthogonal 
basis functions using the singular value decomposition 
such that 



s i(o) =2^Vo)^ u j' 



3i(7T/2) 



(tt/2) 1 



(23) 



(24) 
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FIG. 1. The coefficients Wk used to construct the compos- 
ite detection statistic as a function of basis vector number, 
ordered by their singular value. The solid trace shows the 
coefficients we obtain from (18) assuming a signal amplitude 
A = 20. The dashed trace shows the singular values squared, 
the choice of coefficients given in (10). The dotted trace shows 
"the excess power" choice of coefficients, which are unity for 
all basis vectors. 



This leads to the following expression of the composite 
detection statistic for (22) 

r^^ Wfc (h-u k ) 2 

k 

= ^2 w k Au^Cfc + Bv l ^ /2) (7k + n k ) 



n 2 



(25) 



where rik is a random number drawn from a unit variance 
Gaussian distribution. We have made use of the fact that 
(m • uj) = Sij. 

In order to assess the operating characteristics of (20) 
we simulate several instances of signals and several in- 
stances of noise. We then compare the number of noise 
trials that produce a value of (20) greater than some 
threshold T* with the number of signal, plus noise, trials 
that produce values above the same threshold. This al- 
lows us to parameterize the detection probability versus 
false alarm probability using the value of T* . 



A. Simulations 

In this section we simulate the procedure described 
above. Our goals are the following. First we explore how 
the detection probability of (20) varies of as a function 
of A at a fixed false alarm probability. We verify that 



it peaks when the simulated signal amplitude is equal to 
A. We then compare the performance of (20) with stan- 
dard matched filtering results. We find, to no surprise 
that (20) alone performs worse. However, we also find 
that by using (20) to hierarchically reconstruct the phys- 
ical template SNR the same detection probability can be 
reached for sufficiently low false alarm probability. 

In order to conduct these tests, we apply the singular 
value decomposition to binary neutron star (BNS) wave- 
forms with chirp masses 1.125M Q < M c < 1.24OM and 
component masses 1M@ < m 1 ,m 2 < 3M Q [12]. The 
number of templates required to hexagonally cover this 
range in parameters using a minimal match of 96.8% is 
M = 456, which implies a total number of single-phase 
filters N — 912. These non-spinning waveforms were pro- 
duced to 3.5PN order[16], sampled at 2048 Hz, up to the 
Nyquist frequency of 1024 Hz. The last 10 seconds of 
each waveform, whitened with the initial LIGO ampli- 
tude spectral density, were used to construct the matrix 
of signals S. 

Using this framework we first test how the detection 
probability varies with the choice of A in (20). We simu- 
late ~ 1.8 x 10 7 signals with SNR 7 and evaluate (20) as 
a function of A. Likewise we evaluate (20) for just noise. 
The results of the detection probability at a false-alarm- 
probability of 10~ 3 are shown in Fig. 2. We find that the 
peak of the detection probability occurs when A equals 
the amplitude of the signal, 7. We note that for SNR 7 
signals the low amplitude limit expression (21a) performs 
nearly as well (within a few percent). However, the high 
amplitude limit is considerably worse (by almost a factor 
of two) . SNR 7 was chosen to accentuate the dependency 
of (20) on A. SNR 7 is actually higher than the typical 
SNR threshold one would place in a gravitational-wave 
search. 

In order to test the efficiency of the composite detec- 
tion statistic versus the traditional matched filter ap- 
proach we again generated ~ 1.8 x 10 7 instances of (25) 
for signal and signal plus noise. This time we chose a 
lower, more realistic, signal amplitude of 5. The signals 
had uniform distributions in the template bank and in 
phase angle. We compare this with the standard result 
arising from maximizing the SNR across the bank and 
over phase angle 



[(h- 



3i(0)J 



(h- 



3 i(7T/2) 



(26) 



The result of the two procedures is shown in Fig. 3. 
As expected, the composite detection statistic performs 
worse than explicit reconstruction of the template param- 
eters. However, when the two methods are combined, it 
is possible to reach the same detection probability with 
asymptotically low false alarm probability. 

Fig. 3 contains three curves, showing the detection 
probability Pd versus false alarm probability Pfa when 
there is a signal with amplitude A = 5 in the data. The 
solid-black line is found by choosing the maximum SNR 
across the filter bank as defined in (26). Maximization 
over the filter output is commonly done in gravitational- 
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FIG. 2. Detection probability at a false alarm probability of 
10 -3 as a function of A. We simulated ~ 1.8 x 10 7 signals at 
SNR 7 for 300 values of A in the composite detection statistic 
denned by (20). As we expect the detection probability peaks 
when A = 7, the amplitude of the signal. It is worth noting 
that in this case the low amplitude limit given by (21a) pro- 
vides a similar detection probability. However the large signal 
limit given by (21b) is considerably worse. 
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FIG. 3. Receiver Operator Characteristic curves associated 
with different detection statistics. The solid-black trace shows 
the performance of choosing the maximum likelihood filter 
across the bank. The dashed line shows the performance of 
the composite detection statistic. The dotted line shows the 
performance of choosing max filter across the bank after con- 
ditionally thresholding on the composite detection statistic. 



wave searches. The dashed-black line is the result of 
the composite detection statistic (20) with the choice of 
A = 5. The dotted-black line is the result of first thresh- 
olding on the composite detection statistic and then con- 
ditionally maximizing the SNR over the bank. This pro- 
cedure produces roughly the same detection probability 
for a false alarm probability of 10 -5 but allows one to do 
the full maximization for only 3% of the filtered data. 



B. Use example 

Our simulations indicate that only ~ 3% of the data 
needs to have the physical template parameters recon- 
structed in order to have a similar detection probability 
as the maximum likelihood method at a false alarm prob- 
ability of 10~ J . This section provides an example of what 
this means for a realistic gravitational-wave search. 

Advanced gravitational-wave detectors should be able 
to analyze and locate compact binary sources at the 
moment they merge. Prompt electromagnetic followup 
could confirm a gravitational-wave detection and low- 
latency searches will be critical to maximize the num- 
ber of simultaneously observed signals. However, low- 
latency gravitational-wave searches will be computation- 
ally costly. The reduced filter set proposed by [12] could 
lower the computational cost substantially, helping to en- 
able near real-time searches. Additionally, in this work, 
we have shown that it is possible to reduce the number 
of physical parameter reconstructions by 97% and main- 
tain similar detection efficiencies. If these methods were 
to be used, it would be necessary to understand what the 
result in Fig. 3 implies. 

We now consider how the results of Fig. (3) may ap- 
ply to a low-latency gravitational-wave search and an- 
swer whether or not a false-alarm probability of 10~ 5 is 
a useful operating point. Consider the joint false alarm 
probability for N independent gravitational-wave detec- 
tors in coincidence 



N 



Pi 



fa 



iJV-l 



FA.ii 



(27) 



where Pfa,z is the false alarm probability for the ith de- 
tector and C is the coincidence trials factor. In order to 
understand the double coincidence, limiting false alarm 
rate from the single detector false alarm probability a few 
pieces of information are needed. The first is the num- 
ber of independent trials per second obtained by filtering 
the data. We will take this to be the frequency of the 
minimum point of the noise curve ~ 150 Hz. If one al- 
lows ~ 30 ms to define coincidence this corresponds to 
an additional 5 samples for coincidence trials at 150 Hz. 
Therefore the false alarm rate of double coincidence cor- 
responding to a 10~ 5 false alarm probability in a single 
detector is 150Hz x 5 x 10~ 10 = 7.5 x 10~ 8 Hz = 2.4yr _1 . 
This is well above the false rate that would be required 
for a detection candidate. Therefore we conclude that 



() 



this procedure should not impact the dctectability of near 
threshold signals. 

IV. CONCLUSIONS 

We have presented a study of compact-binary 
gravitational-wave detection that precedes parameter es- 
timation. This could allow more computationally effi- 
cient algorithms to be run in near real time that deter- 
mine whether a signal is present before attempting to 
measure it's parameters. Our study shows that it should 
be possible to reconstruct the physical parameters for 
only of the data while not impacting the sensitiv- 

ity of a compact binary search. 
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